arXiv: 1509.0498Ivl [math.DS] 16 Sep 2015 


A FAMILY OF PERIODIC SOLUTIONS OF THE THREE BODY 

PROBLEM. LIGHT VERSION 


OSCAR PERDOMO 

Abstract. In this paper we describe a 1-dimensional family of initial conditions E that 
provides reduced periodic solution of the three body problem. This family S contains a 
bifurcation point and extend the periodic solution described in [7]. This 1-dimensional 
family is the union of two embedded smooth curves. We will explain how the trajectories of 
the bodies in the solutions coming from one of the embedded curves have two symmetries 
while those coming from the other embedded curve only have one symmetry. The Round 
Taylor Method is a numerical method implemented by the author to keep track of the 
global error and the round-off error. A second version of this paper, same title with the 
“light version” part removed, will include analysis of the error of the solutions using the 
Round Taylor Method. 
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Figure 1. Trajectory of one of the bodies for some of the periodic solutions 
described in this paper 


1. Introduction 

In [7] the author considered the family of solutions of the three body problem described by 
the parameters a and b suggested by the following hgure. Three functions F(t, a, b), R(t, a, b) 
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Figure 2. A family of solutions of the three body problem described by the 
parameters a and b. 

and 0(t, a, b) describe the solution of the three body problem with initial conditions described 
in Figure]^ more precisely we have that 


x{t) = (0,0, F(t, a, 6)) 

y{t) = a, 6) cos(0(t, a, 6)), a, 6) sin(0(t, a, 6)), —F(t, a, 6)) 

z(t) = a, b) cos(0(t, a, 6)), —R{t, a, b) sin(0(t, a, 6)), —Fit, a, b )), 

satisfy the three body problem equations, provided that the following ODE system holds 
true 


( 1 . 1 ) 


F = 


400 


^3 


F, R = 


lOOa^ 

~W~ 


25 


200i? 
^3 ’ 


i?20 = 10a, 


with i?(0, a, b) = 10, F(0, a, b) = 0, 0(0, a, b) = 0, where S = y'R'^ + AF"^ 

where G denotes the partial derivative with respect to t for any function G{t,a,b). We 
are assuming that the gravitational constant is 1. We will be referring to the solution of the 


three body problem described by the differential equation (1.1) as 0(a, b). 

where G denotes the partial derivative with respect to t for any function G{t,a,b). We 
are assuming that the gravitational constant is 1. We will be referring to the solution of the 
three body problem described by the differential equation ( 1.1[ ) as 0(a, b). 

It is not difficult to show that anytime we hnd a point {T,a,b) such that F{T,a,b) = 
R{T,a,b) = 0, then the solution 4’{a,b) is reduced periodic with period 4T. We will call 
these solutions odd/even solutions because from the point of view of f = 0 the function F is 
an odd function with respect to t, but from the point of view of f = T, both functions F and 
R are even. Also, it is not difficult to show that anytime we hnd a point (T, a, b) such that 
F{T, a, b) = R{T, a, b) = 0, then the solution 0(a, b) is reduced periodic with period 2T. We 
will call these solutions odd solutions because from the point of view of f = 0 the function F 
is an odd function. We point out that every odd/even function is also odd due to the fact 
that if F(T, a, b) = R{T, a, b) = 0, then F(2T, a, b) = R{2T, a, b) = 0. 
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In [7j the author proved the existence of a small path a G = {{T, a, h) : T, a, h E R} of 
odd/even solutions. In this paper we extend this path to a path Si of odd/even solutions: 
the embedded curve P 3 BP 2 in Figure]^ This path starts near a periodic solution where the 
body in the center remains motionless and the other two bodies move along an ellipse {b 
near 0). The path ends near a motion with a double coalition (in this case a is near 0). We 
also provide a path S '2 C R^ of odd solutions: the embedded curve that contains the curve 
PiB in Figure This path S '2 starts near a motion with a triple collision (in this case a 
is near 0). This path S '2 intercepts the path S'! in a single point that can be thought as a 
bifurcation point of the set of reduced periodic solutions. The author is not sure if this path 
will continue to be an unbounded curve. 

A second version of this paper will use the Round Taylor Method [7] to compute the error 
and show that the difference between the initial conditions and the values of the solution 
after a period T is small. 

2. Graph oe the path oe initial conditions that provides reduced periodic 

SOLUTIONS. 

This section describes a path C R^ = {(T,a,b) : T,a,b E R} with the property that for 
any (T, a, b) in S, the solution 0(a, b) is reduced periodic with period T, this is, the relative 
position of the three bodies in the solution 0(a, b) repeats every T units of time. Figure is 
an image of the path S. 



Figure 3. Pi has coordinates (5.0063..., 0.109392..., 1.16473 ... ), which 
means that when a = 0.109392..., and b = 1.16473..., then, the solu¬ 
tion 0(a, b) has period T = 5.0063.... The other coordinates are: P2 = 

(12.7012..., 0.0437163..., 3.19541...), P 3 = (9.9472 ..., 4.73605 ..., 0.2) 
and, the coordinate for bifurcation point are B = 
(14.6072 ..., 2.08181..., 3.19493 ...). 

Recall that the distances between the three bodies depend exclusively only on the functions 
R{t) and F{t). More precisely, the distance between the two bodies that go around the z 
axis is 2R{t) and the distance between the body that moves on the z-axis and any of the 
other two bodies is given by ■>/P(t)^ -|- FitY- 

2.1. The solution given by the point P 3 = (9.9472..., 4.73605..., 0.2). This solution 
is very close to the solution when the body in the center stays still and the other two 
bodies move on a perfect circle. Figure shows the graphs of the functions P and R 
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associated with this solution. We can see how the body in the center given by (0,0, F(t)) 
moves very little and the other two bodies given by cos(0(t)), i?(t) sin(0(t)), —F(t)) 
and {—R(t)cos{Q(t)),—R(t)sm{Q(t)),—F(t)) stay near the circle of radius 10 on the x-y 
plane with center at the origin. 



Figure 4. For the point P 3 , the image on the left shows the graph of the 
function F, the image on the center shows the graph of the function R and 
the image on the right shows both functions. 


2.2. The solution given by the point P 2 = (12.7012 ..., 0.0437163 ..., 3.19541...). For 
this motion the body in the center oscillates from (0, 0, —6.01695 ...) to (0, 0, 6.01695 ...). 
This solution is very close to a solution with a double collision. Figure shows the graphs 
of the functions F and R associated with this solution. We can see how after a quarter of a 
period the two bodies that go around the z-axis are very close to each other {R is very small, 
it is near 0.0547972) at this instance the body in the center is at its highest point while the 
other two bodies are both very close to the point (0,0, —6.01695 ...). 



Figure 5. The functions F and R associated with the reduced periodic solu¬ 
tion given by the point P 2 . 


2.3. The solution given by the point Pi = (5.0063..., 0.109392..., 1.16473...). For 
this motion the body in the center oscillates from (0, 0, —1.52497...) to (0, 0,1.52497...). 
This solution is very close to a solution with a triple collision. Figure shows the graphs 
of the functions F and R associated with this solution. We can see how after half of a 
period the three bodies are very close to the origin. F is zero and R is very small, it is near 
0.0369262. 
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Figure 6. The functions F and R associated with the reduced periodic solu¬ 
tion given by the point Pi. 

2.4. The solution given by the point B = (14.6072 ..., 2.08181..., 3.19493 ...). This is 
the bifurcation point of the path of solutions. Figure shows the graphs of the functions F 
and R associated with this solution. 



Figure 7. The functions F and R associated with the reduced periodic solu¬ 
tion given by the point B 

Let us call Si, the part of the path contained in S given by the smooth curve that connects 
P 2 , B and P 3 , see Figure Likewise, let us call S 2 , the part of the path contained in S 
given by the smooth curve that connects Pi, B and also contains the part of S that has an 
arrow at the end. All the points in the path Si are odd/even solutions; all of them share the 
symmetry given by the functions shown in Figure More precisely: With respect to the 
origin, the function P is odd and the function R is even and with respect to t = both 
functions P and R are even. The solutions associated with points in the curve S 2 different 
from the point B do not share these two symmetries. With respect to t = 0, the function P 
is odd and the function R is even. There is not symmetry with respect to t = Figure]^ 
shows the function P and R for the solution associated with two points on 5 * 2 . 

3. Progedure to obtain the points on the gurve S and the bieurgation point 

3.1. Original ODE. The functions F{t,a,b), R{t,a,b) and Q{t,a,b) satisfy the following 
ODE 
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Figure 8. The image on the left shows the functions F and R associated with 
a point on the curve S 2 between Pi and B. The image on the right shows these 
two functions for a point on the curve S 2 located after the point B. We can 
see how the symmetry with respect to the vertical line ^ f that is present 
for solutions coming from points in Si is not present in these solutions. 


Xi = X-i 

X2 = Xa 

400x1 
(4xi + X 2 ) 

. _ lOOa^ 25 200x2 

X 2 X 2 (4xf + X 2 ) 
10a 


with initial conditions 

Xi(0) = 0 X2(0) = 10 X3(0) = b X4(0) = 0 X5(0) = 0 

This ODE will be referred to as the original ODE. We have that Xiit) provides P(f, a, 6), 
X 2 (f) provides R(t,a,b), x^it) provides F(t,a,b), x^it) provides R(t,a,b) and x^it) provides 
0(f, a, b) 

3.2. Extended ODE. The partial derivative of the functions F{t,a,b), R{t,a,b) and 
0(f, a, b) satisfy the following ODE 


Xi 

X 2 

X 3 

X 4 

X 5 

xe 

X 7 

Xs 


X 3 

X 4 

400xi 

(4x1^ + X 2 ^) 

200x2 _ 25 lOOa^ 

(4x1^ + X 2 ^) X 2 ^ X 2 ^ 

10a 

X2‘^ 

Xs 

X 9 

3200x6Xi^ 1200x2X7X1 400x2^X6 

(4xi2 + X 2 ^) 5/2 + (4xi2 + x22) ^2 ~ (4a:,2 + 
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3000^X7 200 a 2400 xiX 2 Xq 50 x 7 600a;2^a;7 200x7 

X2^ X2^ (4X1^ + X2^) X2^ (4X1^ + X2^) (4X1^ + X2^) 

10 20 ax 7 

X2^ X2^ 

Ti3 

Xi4 

3200xiiXi^ 1200 x 2 X 12 X 1 400x2^X11 

(4x1^ + X 2 ^) (4xi^ + X 2 ^) (4X1^ + X 2 ^) 

300a^xi2 2400x1X2X11 50 xi2 600x2^xi2 200xi2 

X 2 "‘ (4x1^ + X 2 ^) ^ X 2 ^ (4x1^ + X 2 ^) (4X1^ + X 2 ^) 

20 axi 2 


with initial conditions 

xi(0) = 0 
X2(0) = 10 
xsjo) = b 

X 4 ( 0 ) = 0 

X5(0) = 0 

X6(0) = 0 

X7(0) = 0 

X8(0) = 0 

X9(0) = 0 

xio(O) = 0 

Xii(O) = 0 
0 : 12 ( 0 ) = 0 

Xi3(0) = 1 

0 : 14 ( 0 ) = 0 

Ti5(0) = 0 

This ODE will be referred to as the extendded ODE. We have that Xi(t) provides F(t, a, 6), 
X2(t) provides R(t,a,b), x^it) provides F(t,a,b), X4(t) provides R(t,a,b), x^{t) provides 
Q{t,a,b), Xfiit) provides ^{t,a,b), X7(t) provides ^{t,a,b), xs{t) provides ^^{t,a,b), 
xg{t) provides ^^{t,a,b), Xio(t) provides ^{t,a,b), Xii(t) provides ^{t,a,b), Xi2(t) pro¬ 
vides ^{t,a,b), Xis{t) provides ^^{t,a,b), Xi4(t) provides ^j^{t,a,b) and Xi5(t) provides 
^{t,a,b). 

These equations easily follows as an application of the chain rule. To exemplify this 
procedure let us compute X8(^); 

Xs = Xia = (-400X1 (4xi^ -1- X2^) 

= —400x6 (4xi^ -I- X 2 ^) -I- 1200xi (4xi^ 4- X 2 ^) (4xiX6 + X 2 X 7 ) 

3200 x 0 X 1 ^ 1200x2X7X1 400 x 2 ^X 6 

(4xi2 + X22) 5/2 + (4 xi2 + X 22 ) 5/2 “ (4xi2 + X 22 ) 5/2 


Xg = 

a^io = 

Xii = 

T 12 = 
T13 = 

Xi4 = 

3^15 = 
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3.3. Getting the points on Si. Ideally we want each point {t,a,b) G Si to satisfy the 
system of eqnations {F(|,a, 6 ) = R{^,a,b) = 0}. From the paper [7] we know that there 
exist a point Pq = (t,a, b) near 


Pq — {to, 0 , 0 , bo) — 


/13366894627923 43170475352787 1490359743\ 


V 5000000000000 ’ 10000000000000’1000000000 J 


snch that F{Po) = R{Po) = 0 . If we assnme for a moment that we completely know the 
fnnctions F and R in the whole R^, then the way to hnd the curve Si would be simply, hrst, 
hnd the vector held X = vi x V 2 , where Vi = VF = {F, Fa, Ff,) and V 2 = Vi? = {R, Ra, Rb) 
and second, hnd the integral curve of the vector held X that goes through Pq. Recall that 
the hrst entry of points in Si is 4 times the hrst entry of points in this integral curve. Notice 


that using the notation from section 3.2 we have that 


VF = 


400x1 (t) 


(4xi(f)2 + X2{t)‘^) 


, Xi3 


and 


Vi? 


25 


8 x2 (t)^ 


(4xi(t)2+X2(t)2)3/2 


X2{t) + 4a^ 


X 2 {ty 


, Xo,Xii 


Even though we do not know the vector held VF and VF everywhere, we can still use the 
idea of integrating the vector held X using the Euler method. The next subsection explains 
the algorithm that we are using in this paper to get the points in Si. 


3.3.1. Continuation algorithm. For the Algorithm we select three small numbers ei and 62 
and €3 that we use as tolerance for the error. We want 0 < ei < € 2 . Recall that we want to 
hnd solutions of the system {F = F = 0}. We will be collecting the solution of this system 
of equations in a set called TS. To start with, we make, TS = {Fq}, the set which only 
element is the numerical solution that we know. We will use Pq to start the algorithm. 


( 1 ) 

( 2 ) 

( 3 ) 


( 4 ) 

( 5 ) 

( 6 ) 
( 7 ) 


Consider a point Qq = {i,a,b) such that |F((5o)| < ci and |F(Qo)| < ^i- 
Make Y = Qq and select a positive real number h and a positive integer k. 

For i = 1 to i = k Do 

• Find the values of X by numerically solving the extended ode (see section 3.2) 
using the values of t, a and b given by the entries of Y. This is, we integrate for 
t units of time where t is the hrst entry of Y and we use the values of a and b 
given by the second and third entry of Y to provide the initial conditions of the 
ode. 

• Make V = V + ^ X and g* = V. 


If for all i = l...k, |F(gj)| < 62 and |F(gj)| < €2 then go to the next numeral, 
otherwise go back to numeral ( 2 ) and select a smaller h and/or a smaller k. 

Find Qi such that |F(Qi)| < ei and |F(Qi)| < ei and such that it is near in the 
sense that \qk — Qi\ < £3 
Make TS = TSC {Qi} U 

Go back to numeral ( 1 ) to start the process over by making Qq = Qi. 


Definition 3.1. We will call pillar points all the points Qi in the previous algorithm. 
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In our case the algorithm stopped or became difficult to carry near the point P 3 and P 2 
because the ode was approaching a singularity near P 2 (recall that P 2 is close to a collision) 
and on the other hand the vector held X is close to the zero vector near P 3 . 

As an example of the method, if we start with Qq = Pq and we use ei = 10 “®, €2 = = 

0.00005, h = 0.001 and k = 200, then, 
gaoo = (2.7219062659312807,4.212655007080421,1.6538674269975053) 
and we can use Qi = (2.72191575576588,4.212633490447383,1.6538497779324066). 


3.4. Getting the points on Si. Ideally we want each point {t,a,b) G S 2 to satisfy the 
system of equations {P(|,a, 6 ) = R{^,a,b) = 0}. If we assume for a moment that we 
completely know the functions P and R in the whole R^, then the way to hnd the curve S 2 
would be simply, hrst, hnd the vector held Z = wi x W 2 , where Wi = VP = (P, Fa, Ft) and 
W 2 = VP = (P, Pa, Rb) and second, hnd the integral curve of the vector held Z. Recall that 
the hrst entry of points in S 2 is twice the hrst entry of points in this integral curve. Notice 
that using the notation from section 3.2 we have that VP = {x^,xq,xii). The points in S 2 


are found using the algorithm described in section 3.3.1 


3.5. Getting the bifurcation point. Recall that the solutions associated with points on 
Si are called odd/even solutions and those associated with points in S 2 are called even 
solutions. Having in mind the symmetries of the solutions, it is not a surprise that for points 
on the integral curve 7 of the vector held X that passes to the point Pq, the vector held Z is 
parallel to the vector held X. The bifurcation point occurs because along 7 , there is a point 
where the vector held Z vanishes. In order to numerically compute the bifurcation point, 
we take the minimum of the set {|Z(|,a, 6 )| : {t,a,b) G S'!}. This minimum occurs at the 
point B = (14.607249047056753,2.081806260749908,3.194934273913219) G ^ 1 . 


4. Periodic Solutions 

Recall that a reduced periodic solution b) with period T is periodic if is a 

rational number. It is not difficult to see the values of the function 0 on the points on the 
S. In this section we will be selecting some values for 0 along the path S that provides 
periodic solutions. 


4.1. Periodic solution displayed on Figure Every entry of the matrix in Figure [T] 
(located at the begging of the paper) shows the trajectory of one of the three bodies that 
moves periodically solving the three-body problem. The following table provides the initial 
conditions of these 9 periodic solutions. The hrst entry T represent the period of the reduced 
periodic solution. This is, the positions and velocities after T units of time agree with the 
initial positions and velocities up to a rotation of 9 radians. This rotation 9 is given in the 
second column of the table. Every color in these images show the trajectory after T units 
of time. These points are part of the path S' 2 . These periodic solutions only one symmetry. 
The periodic solution given by the last row in this table is very close to a triple collision. 
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Coordinates: (T, a, b) 

Q{T,a,b) 

Image 

location 

(7.464725167070125, 1.332130235206886, 2.39131806234605) 

147r 

9 

(1,1) 

(6.692611549615348, 1.1189865713077587, 2.178542324667063) 

iItt 

7 

(1.2) 

(6.156645197408822, 0.9203809141218081, 1.9818897189116862) 

5 

(1.3) 

(5.758992220584509, 0.72863210825991, 1.791870916952545) 

l^TT 

8 

(2.1) 

(5.454531817007846, 0.5397530375918577, 1.6029431700479464) 

bTT 

3 

(2.2) 

(5.32949878269853, 0.4457220498597047, 1.5077112808284203) 

VItt 

7 

(2.3) 

(5.220379172126002, 0.35201283725645105, 1.411839280497764) 

/TT 

4 

(3.1) 

(5.1264630045948305, 0.25902392732426605, 1.3159094056731593) 

IbTT 

9 

(3.2) 

(5.096931182326409, 0.22679159236240487, 1.2822928944673215) 

977 

_5_ 

(3.3) 


4.2. Periodic solution displayed on Figure Every entry of the matrix in Figure 
1 ^ (located at the end of the paper) shows the trajectory of one of the three bodies that 
moves periodically solving the three-body problem. The following table provides the initial 
conditions of these 9 periodic solutions. The hrst entry T represent the period of the reduced 
periodic solution. This is, the positions and velocities after T units of time agree with the 
initial positions and velocities up to a rotation of 6 radians. This rotation 6 is given in the 
second column of the table. Every color in these images show the trajectory after T units of 
time. These points are part of the path Si. These periodic solutions are odd/even solutions 
and therefore they have two symmetries. 


Coordinates: (T, a, b) 

Q{T,a,b) 

Image 

location 

(10.694782129146047, 4.3162773916465715, 1.4916623030663265) 

17 TV 

(1,1) 

(10.900907564917922, 4.205530359018152, 1.6642665366638942) 

i^TT 

7 

(1.2) 

(11.251546754899636, 4.020933846016405, 1.9098147447184282) 

IbTT 

6 

(1.3) 

(11.546587626864484, 3.8684643701482133, 2.0840362602898437) 

IItt 

5 

(2.1) 

(12.005866456188725, 3.6348152391561275, 2.314439960284758) 

9^ 

4 

(2.2) 

(12.482298481771874, 3.3941653076488447, 2.5161535540867117) 

ibvr 

7 

(2.3) 

(12.805950195048542, 3.229502929775678, 2.6373059506494907) 

Ytt 

3 

(3.1) 

(13.037980888481863, 3.10968118474065, 2.717938296845948) 

iZTT 

5 

(3.2) 

(13.211458502729283‘, 3.01851140163578', 2.775360523910528') 

1(7V 

_1_ 

(3.3) 


4.3. Periodic solution displayed on Figure |10|, Every entry of the matrix in Figure 


10 (located at the end of the paper) shows the trajectory of one of the three bodies that 
moves periodically solving the three-body problem. The following table provides the initial 
conditions of these 9 periodic solutions. The hrst entry T represent the period of the reduced 
periodic solution. This is, the positions and velocities after T units of time agree with the 
initial positions and velocities up to a rotation of 6 radians. This rotation 6 is given in the 
second column of the table. Every color in these images show the trajectory after T units of 
time. These points are part of the path Si. These periodic solutions are odd/even solutions 
and therefore they have two symmetries. 
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Coordinates: (T, a, h) 

0(T,a,6) 

Image 

location 

(13.451787406253272, 2.8889144618608733, 2.8514754095228265) 

IItt 

6 

(1,1) 

(13.60916416956247, 2.8011806024112604, 2.8994928664781003) 

IcStt 

7 

(1.2) 

(13.719567754906889, 2.737827736709047, 2.9324647582021814) 

iOTT 

8 

(1.3) 

(13.801004797570164, 2.689923890770301, 2.9564685196206826) 

i/TT 

9 

(2.1) 

(14.320649658996734, 2.344448198979306, 3.106748248260848) 

271 

(2.2) 

(14.657003574778068, 2.0202122629047206, 3.212388918731313) 

IYtt 

l^TT 

(2.3) 

(14.686554119081652, 1.9783743950165875, 3.2235432541405697) 

7 

(3.1) 

(14.719531067694582, 1.9241180805387452, 3.2371600512994565) 

i^TT 

6 

(3.2) 

(14.754026030069339, 1.8509347720107878‘, 3.254002265714927) 

ilTT 

_5_ 

(3.3) 


4.4. Periodic solution displayed on Figure |11|, Every entry of the matrix in Figure 


11 (located at the end of the paper) shows the trajectory of one of the three bodies that 
moves periodically solving the three-body problem. The following table provides the initial 
conditions of these 9 periodic solutions. The hrst entry T represent the period of the reduced 
periodic solution. This is, the positions and velocities after T units of time agree with the 
initial positions and velocities up to a rotation of 6 radians. This rotation 6 is given in the 
second column of the table. Every color in these images show the trajectory after T units of 
time. These point are part of the path Si. These periodic solutions are odd/even solutions 
and therefore they have two symmetries. 


Coordinates: (T, a, h) 

0(T,a,6) 

Image 

location 

(14.782277611145048, 1.7467887078517095, 3.274928737204819) 

yTT 

4 

(1.1) 

(14.786834114569135, 1.6762115129006632, 3.287061516889608) 

ItlTT 

7 

(1.2) 

(14.774959957278682, 1.5866395627104857, 3.300053407198923) 

Ytt 

3 

(1.3) 

(14.7286222583901, 1.4691728012133716, 3.3129725743424996) 

IZTT 

5 

(2.1) 

(14.699688133885457, 1.4214348721248662, 3.316877441047812) 

IYtt 

7 

(2.2) 

(14.60791563192398, 1.3083245035642173, 3.3230110738911285) 

57r 

2 

(2.3) 

(14.4942511626078, 1.2033604318255704, 3.324782465355579) 

7 

(3.1) 

(14.444161869664121, 1.163467447374621, 3.3244726259938155) 

137r 

5 

(3.2) 

(14.319901500300364, 1.0746785092812807, 3.321867812104126) 

^TT 

_3_ 

(3.3) 


4.5. Periodic solution displayed on Figure \12\ Every entry of the matrix in Figure [T^ 
(located at the begging of the paper) shows the trajectory of one of the three bodies that 
moves periodically solving the three-body problem. The following table provides the initial 
conditions of these 9 periodic solutions. The hrst entry T represent the period of the reduced 
periodic solution. This is, the positions and velocities after T units of time agree with the 
initial positions and velocities up to a rotation of 9 radians. This rotation 9 is given in the 
second column of the table. Every color in these images show the trajectory after T units of 
time. These point are part of the path Si. These periodic solutions are odd/even solutions 
and therefore they have two symmetries. The periodic solution given by the last row in this 
table is very close to a double collision. 
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Coordinates: (T, a, h) 

0(T,a,6) 

Image 

location 

(14.155328130457113, 0.9714850334241202, 3.3156069432644326) 

IItt 

4 

(1,1) 

(14.054387879818348, 0.9133674393862702, 3.310617589529648) 

147r 

(1.2) 

(13.986985590840762, 0.8760917450402825, 3.306890869738851) 

iVtt 

6 

(1.3) 

(13.659851940872658, 0.70551977441241, 3.285155344523193) 

37r 

(2.1) 

(13.315322141213205, 0.5298843477532729, 3.2569218042195365) 

lt)7r 

.5 

(2.2) 

(13.240214962690887, 0.4900089555434633, 3.250127637396059) 

loTT 

4 

(2.3) 

(13.125795380869723, 0.4265448292958072, 3.239357536186598) 

TDtf 

3 

(3.1) 

(13.0441592363555, 0.3781732129466524, 3.23136706052875) 

l/TT 

5 

(3.2) 

(12.938357527438528, 0.30903261362830825, 3.2206342443812517) 

Ytt 

4 

(3.3) 
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Figure 9. Trajectory of one of the bodies for some of the periodic solutions 
described in this paper. 
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Figure 10. Trajectory of one of the bodies for some of the periodic solutions 
described in this paper. 
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Figure 11. Trajectory of one of the bodies for some of the periodic solutions 
described in this paper. 
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Figure 12. Trajectory of one of the bodies for some of the periodic solutions 
described in this paper. 



















